YITP- 12-66 



(N 

o 

(N 

> 

O 

O: 
ov 



o 
cr 



Impact of the second order self-forces on the dephasing of the gravitational waves from 

quasi-circular extreme mass-ratio inspirals 

Soichiro Isoyama 1 !^ Ryuichi Fujita 2 Norichika Sago 3 Hideyuki Tagoshi 4 J§ and Takahiro Tanaka^ 

1 Yukawa Institute for Theoretical Physics, Kyoto university, Kyoto, 606-8502, Japan 
2 Departament de Fi'sica, Universitat de les Illes Balears, Palma de Mallorca, E-07122 Spain 
3 Faculty of Arts and Science, Kyushu University, Fukuoka 819-0395, Japan 
4 Department of Earth and Space Science, Graduate School of Science, Osaka univeristy, Osaka, 560-0043, Japan 

(Dated: October 10, 2012) 

The accurate calculation of long-term phase evolution of gravitational wave (GW) forms from extreme (inter- 
mediate) mass ratio inspirals (E(I)MRIs) is an inevitable step to extract information from this system. Achieving 
this goal, it is believed that we need to understand the gravitational self-forces. However, it is not quntatively 
demonstrated that the second order self-forces are necessary for this purpose. In this paper we revisit the prob- 
lem to estimate the order of magnitude of the dephasing caused by the second order self-forces on a small body 
in a quasi-circular orbit around a Kerr black hole, based on the knowledge of the post-Newtonian (PN) approx- 
imation and invoking the energy balance argument. In particular, we focus on the averaged dissipative part of 
the self-force, since it gives the leading order contribution among the various components of them. To avoid 
the possibility that the energy flux of GWs becomes negative, we propose a new simple resummation called 
exponential resummation, which assures the positivity of the energy flux. In order to estimate the magnitude 
of the yet unknown second order self-forces, here we point out the scaling property in the absolute value of the 
PN coefficients of the energy flux. Using these new tools, we evaluate the expected magnitude of dephasing. 
Our analysis indicates that the dephasing due to the second order self-forces for quasi-circular E(I)MRIs may 
be well captured by the 3PN energy flux, once we obtain all the spin dependent terms, except for the case with 
an extremely large spin of the central Kerr black. 
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I. INTRODUCTION 

Extreme mass ratio inspirals (EMRIs) and intermediate mass ratio inspirals (IMRIs), in which a stellar mass or several tens 
solar mass compact object inspirals into a more massive central black hole, have attracted much interest not only as a promising 
source of the gravitational waves (GWs) for future space born GW detectors, but also as a unique clean probe of the spacetime 
region of strong gravity. To achieve the test of general relativity using GWs from E(I)MRIs, we need to predict sufficiently 
accurate waveforms. This requirement motivates us to model the E(I)MRIs as the motion of a small body in a given backgroud 
spacetime with gravitational back reaction. This back reaction is treated as gravitational self-forces ijj-jg] and its higher order 
extension with respect to the mass ratio has attracted much interest in recent years 0-Q1]]. See following review articles, e.g. 
, Il2l[l3ll and references therein for more details. 

In fact, following the scaling argument (See e.g. fjilfislp , the phase of GWs from a particle whose orbital frequency sweeps 
a few orders of magnitude before the plunge can be expanded as 
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where <f>(°' and are 0(1) quantities independent of fx, which is the mass of a small particle, and Mbh and a are the mass 
and the spin parameter of the Kerr black hole, respectively. 1 On the one hand, in Eq. ([1]), we only need the self-forces 
up to the first order time-averaged dissipative part. In addition, it has been long known that can be computed with well 
established balance argument, which relates the first order time-averaged dissipative part of self-forces to the energy and angular 
momentum fluxes associated to the global killing vectors on background Kerr black hole lfl9ll . Though the Carter constant is 
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not associated to any Killing vector and simple balance argument is not applicable, the way to computing long time averaged 
evolution of Carter constant has been already well established 1I20I - I23I1 . On the other hands, <5>^\ which is referred to as the 
dephasing relative to <E>(°), depends on two different components of the self-forces: the first order conservative part, and the 
averaged second order dissipative part. 2 The first sub-leading term still can be important since potentially it may give the 
correction significantly greater than unity to the phase 124142611 . 

To obtain accurate waveform, there have been many works on self-forces. As for the first order conservative part, thanks to 
the recent massive development, we are now in part ready for practical computation with numerical implementation I2~7l - t3lll . 
Particularly, in the case of Schwarzschild background the correction to the orbital frequencies l32ll and lO have been 
already studied extensively. Even in the case of Kerr background the preliminary results of the self-forces in quasi-circular 
orbits have been reported 113411 . 

By contrast, the averaged second order dissipative part of self-forces has been so far studied only at the formal level in the 
context of black hole perturbation l7l 4Tlll . It will require much more effort to establish the method to compute the second order 
dissipative part, especially in the case of Kerr background. Under such circumstances, a typical strategy to evaluate is 
to make use of the standard post-Newtonian (PN) approximation, in which we assume slow motion of a satelite and its weak 
gravitational field. Based on the PN approximation, Huerta and Gair IB5I1 evaluated the size of dephasing caused by the first order 
conservative self-forces and the averaged second order self-forces, pickin g up representative EMRls in quasi-circular orbits on 
a Kerr black hole. The same dephasing was also discussed by Yunes et al. [36] using the effective one body formalism again for 
representative E(I)MRIs. 

A naive expectation is that the PN approximation is not suitable for modelling the waveforms of E(I)MRIs, especially for a 
Kerr black hole with large spin. Typical E(I)MRIs in circular orbits spend the last few years of inspiral in the vicinity of the 
inner most stable circular orbit (ISCO). Since the ISCO radius reaches the event horizon as the spin of the Kerr black hole is 
increased to the extremal limit, the motion of the body becomes highly relativistic, exceeding the validity range of the standard 
PN approximation l37ll . However, the things are not so trivial. The time spent near the ISCO becomes longer as we increase the 
mass ratio Mbh/m> but in that case the mass of the large central black hole also becomes larger. Then, the total cycles of GWs 
become smaller for a given observation time. As a result, the correction due to higher order self-forces might be suppressed 
below the observational threshold, despite the loss of the accuracy of the PN approximation. Therefore, it is not so obvious if 
there are E(I)MRIs such that really require the notion of the second order self-forces. 

The previous analyses mentioned above JH, are focusing on the corrections coming from the self-forces at the currently 
available PN order and are limited to representative E(I)MRIs. To get an insight into whether the second order self-forces based 
on the black hole perturbation are really necessary to calculate the waveforms of quasi-circular E(I)MRIs or not, therefore, it 
would be useful to give an estimate of dephasing coming from the averaged second order dissipative self-forces, focusing on the 
yet unknown higher PN terms and surveying the whole parameter region of E(I)MRIs 

What we discuss in this paper is the adiabatic evolution of E(I)MRIs in quasi-circular orbits on a Kerr spacetime. Here the 
adiabatic evolution means an approximation in which the evolution of the orbital frequency is determined by the energy balance 
argument, i.e., the rate of change of the total energy of the binary is equated to the energy flux emitted to the infinity. 

To evaluate the order of magnitude of the yet unknown higher PN corrections, we need to rely on some extrapolation. For this 
purpose, we first introduce a simple new resummation of the energy flux, which we call the "exponential resummation". When 
the spin of a Kerr black hole is large enough, the PN energy flux in the Taylor form can be negative outside the ISCO radius for 
some PN orders j38ll . If this happens, the estimated total phase before the plunge diverges and the extrapolation to the higher 
PN order will not make sense. Our simple "exponential resummation" is the one that ensures the positivity of the energy flux. 

As a PN input for the corrections at the next leading order in the mass ratio, the best one available so far is the 3.5PN energy 
flux of GWs II39L [4011 with linear spin dependent terms up to 3PN order, which has been recently derived by Blanchet et al. 14111 . 
To estimate the possible magnitude of the yet unknown higher order PN terms, we focus on a scaling property among the PN 
coefficients in the energy flux. Using the 8PN energy flux in the test particle limit l42l l43ll . we will show that the absolute 
values of the coefficients scale roughly as required from the convergence of the PN series up to the light ring radius. (This point 
was also discussed by Nakano et al. independently ll44ll .) Since this scaling behavior is related to the PN convergence, despite 
the lack of the higher PN terms in the energy flux, we conjecture that the same scaling property will hold for the terms higher 
order in the mass ratio. Under this assumption, we estimate the order of magnitude of the unknown portion of the energy flux 
coming from higher PN terms at the next leading order in the mass ratio via "extrapolation". Gathering the tools stated above, 
we will investigate the impact on the dephasing from the averaged dissipative second-order self-forces for various E(I)MRIs 
systematically. This is the main goal of this paper. 



2 Here the "dissipative" part means the self-forces that cause the time variation of the constants of motion, such as the energy, the angular momentum around 
the axis of symmetry and the Carter constant. The "conservative" part means the one giving the correction to the relation between the orbital frequencies and 
the constants of motion fT3 . l2fj|| . The meaning of "time-averaged" is averaging over a sufficient long period compared to the time scale for the evolution of 
the phase difference between the oscillations in the radial and the zenith angle directions. 
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The remainder of this manuscript is organized as follows. In Sec. [II] we briefly review how the accumulated phase of GWs 
from adiabatic inspiral is calculated, based on the balance argument. In Sec. [TTTJ to cure the negative energy flux of GWs that 
appearers in the truncated PN Taylor series expansion, we propose the exponential resummation that ensures the positivity of 
the energy flux. Sec. [IV] is dedicated to the study of the scaling property in the coefficients of the energy flux in the test particle 
limit that becomes manifest owing to the brand new 8PN energy flux |43J|. In Sec.[V] using the exponential resummation and the 
scaling property, we will estimate the dephasing coming from the yet unknown part of the second order self-forces for various 
E(I)MRIs. We find that the unkown non-linear spin dependent terms at the lower PN order in the energy flux dominante the 
unknown dephasing. We summarize our results and conclude in Sec.lVIl 

In this manuscript, we use geometrical units G = c = 1 and the sign convention of the metric is (—,+,+,+). The coordinates 
(t, r, 9, (f) denote the Boyer-Lindquist coordinates of the Kerr black hole l45ll . We frequently use the dimensionless spin defined 
by q := a/MsH, and the symmetric mass ratio defined by v := \iMbhI (Mbh + m) 2 - 



II. THE ACCUMULATED PHASE OF THE GRAVITATIONAL WAVE FROM AN INSPIRALING BINARY 

We consider a binary composed of a small satellite body with the rest mass /i in a quasi-circular orbit around a Kerr black 
hole with the mass Mbb and the spin parameter a. We neglect the effect of the spin of the satellite, which is negligiblly small 
for E(I)MRIs though plays no negligible role in some situations ll46ll47tl . For a binary in a quasi-circular orbit, the accumulated 
phase of GWs is calculated as 

r x ° x 3/2 e'(x) 

$ := -2 / dx—t^ , (2) 

J xlsco M E(x) 

where E is the binding energy of the binary, E := dE/dt is the energy loss rate, and the prime denotes the differentiation with 
respect to x, which is a dimensionless orbital frequency defined by x := (Mf2) 2 / 3 with the orbital frequency fl and the total 
mass M := Mbh + Here E is supposed to depend only on x, neglecting the effect of variation of mass due to the energy 
absorbed by the black hole. Throughout the paper we totally neglect the effect of this time dependent mass variation. We choose 
the lower bound of the integral in Eq. (O, xisco to the value of x at the inner most stable circular orbit (ISCO) and the upper 
bound xq to the one determined by the condition coming from the finite observation time. Throughout this paper, we adopt 
£obs = 1 year as the observation time. 

Our primary interest in this paper is the order of magnitude of the dephasing coming from the corrections at the next leading 
order in the mass ratio. To evaluate it, we calculate the difference between $ with and without the higher order corrections. 
Checking the detectability of this difference by observations, the initial and the final frequencies should be kept unchanged when 
we evaluate it. Here, we fix both xq and £isco to the values determined by the test particle limit. Namely, 114811 



xisco := ( i?isco + 1 ) ' ftsco := 3 + Z 2 T y/@ - Zi)(3 + Z x + 2Z 2 ) , 



-1/6 



Zi := 1 + (1 - <Z 2 ) 1/3 [(1 + q) 1/3 + (1 - q) 1/3 } , Z 2 := y/3q 2 + Z\ , (3) 

where the upper (lower) sign in the second equation is chosen for the co (counter)-rotating case with q > (<)0 and R := t/Mbb 
represents the dimensionless radius, xq is determined by 

1 (year) = - / dx . , / ; , (4) 

where the superscript I°l means the leading order contribution in the mass ratio, i.e., the test particle limit, and 

m _ R^-2R^ + g 1 p _ f l-qx^ \ V3 



'#3/4(^3/2 _ 3J jl/2 + 23)1/2 ' ■ Af(i? 3 / 2 + g)' V 

As for the energy, we here adopted the binding energy of the binary instead of the energy directly related to the four momentum 
of the particle. It would be more natural to consider the latter in the context of the self-force calculation based on the black hole 
perturbation, but it is not gauge invariant. Hence, the former is more suitable for comparison with the PN calculation. These two 
energies are different at the next leading order in the mass ratio due to the presence of the gravitational field energy, but there 
must be one-to-one correspondence between them, once the gauge is completely fixed. 

The energy balance argument tells us that the averaged loss of the total energy should be equal to the averaged total energy 
flux to the future null infinity because of the energy conservation (see e.g. Ref. IU9H49ll ). Hence, as long as we define the binding 
energy appropriately, 
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holds, in the sense averaged over a sufficiently long time and neglecting the horizon absorption flux, where £ is the energy flux 
emitted to the null infinity. Hence, we evaluate the phase by 

In the above equation, as we are interested in the dissipative corrections, we fixed E' to the expression in the test particle limit. 
We expand £ as £ = +vL^ +0(t^). To compute we can use the Teukolsky formalism ||50|| and invoke the numerical 
code developed by Fujita and Tagoshi l5l"l l52tl . 

The finite mass corrections at the next leading order, are in part provided by the standard PN calculations. In the PN 
formalism, the energy flux of GWs emitted to the infinity from a quasi-circular binary is obtained up to 3.5PN order for the 
spin independent terms lf39ll40ll and up to 3PN order for the terms linear in spin when the spin vectors are parallel to the orbital 
axis |41j. In the present notation, truncated at 0{v), it is given by 

/ 1247 35 \ , /2 / 11 17 \ 2 / 44711 9271 

1 + x [ v + x i,z 47r q H qv } + x 2 [ 1 V 

336 12 J \ 4 q 4 y J \ 9072 504 



e,a f 8191 59 f 583 3749 1 \ 

+x b/2 7T q + v { 7r H q \ 

\ 672 16 y \ 24 144 q J J 



,( 6643739519 16 , 1712 856, . , 65 f 134543 41 2 33 

+X { 69854400 + ~ T05~ 7E ~ 105 lo ^ 16x ^ _ ~Q 1Tq + ( 7775" + 48 ^ + 

7/9 / 16285 214745 \ , , , 7/9 " 

+x 7 ' 2 1 v it + O x 4 , v 2 , q 2 , qx 7/2 ) 

\ 504 1728 J K ' ' y ' 



(8) 



where 7e = 0.57721 • • • denotes Euler's constant. We will refer to the above expression for the energy flux truncated at the nth 
PN order 3 as "the nPN Taylor flux". We denote it by £^ PN , and expand it like £^ PN := £^pn + v^nVN + " " as before. We 
also frequently use "the normalized nPN flux" defined by 



£nPN<>,<?) := (^y^ 5 ) £ nPN(>>9)- (9) 

From now on, we denote the PN flux solely composed of the known part as "the known nPN Taylor flux" and distinguishing it 
from the full PN flux by associating " ~ " like £^ PN . We also introduce the notation to denote the residual terms higher order in 
PN expansion and in spin by 

r T[i] ._ r [i] _ f T[i] nm 
A '>riPN - — ^full ^nPN- { Lv > 

~ Till 

Then, L >3 5PN represents the yet unknown energy flux to be determined from the computation of the second order dissipative 
self-forces. Further we define <C^pjj and £>^ PN by the terms in that correspond to £^pjj and £>'^ PN , respectively. To be 
precise, we define £ raP jj by the sum of the spin independent terms up to 3.5PN order and the terms linear in the black hole spin 
up to 3PN order. £ > , lPN is the remaining terms at the leading order in the mass ratio, which also includes the non-linear spin 
dependent terms in all PN orders. 4 . 



III. A MANIFESTLY POSITIVE DEFINITE ENERGY FLUX: EXPONENTIAL RESUMMATION 



In this section we propose a new simple resummation scheme, which can be easily applied if we just know the nPN Taylor 
flux. 



A. Limitation of the Taylor flux revisited 

It is easy to imagine that the PN Taylor flux is not so accurate when the orbital radius becomes small. Indeed, Zhang et al. i37ll 
pointed out that the nPN Taylor flux in the test particle limit rapidly loses the accuracy around the ISCO radius. However, it is a 



3 We refer to a term of O (x n ) relative to the leading order as those of nPN order. 

4 Here we do not distinguish the logarithmic term in 3PN Taylor flux in the test particle limit from the other 3PN terms. We confirmed that the results in the 
present paper do not change much even if we exclude this term from •C^p^ 
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FIG. 1: The normalized nPN Taylor flux in the test particle limit, L^pl,(x, q), for q = 0.9 up to £isco(0.9) := 0.78014 .... The 
horizontal axis is the dimensionless frequency x of a body. The lable "Exact" means the exact numerical energy flux, calculated by Fujita and 
Tagoshi ifsTh . 



different issue whether or not the Taylor flux is accurate enough for our present purpose because what we are interested in here 
is whether or not the dephasing in GW waveforms is measurebale. 

As we want to evaluate the correction to $ in Eq. (O caused by the yet unknown part of the flux £>„p N , we need some 
extrapolation method. For this purpose, later it becomes necessary to evaluate something like the phase for the flux truncated 
at nPN order, $[£^p ^]. At this point we find that the nPN Taylor flux is problematic. As mentioned earlier, the energy flux 
truncated at some PN orders become negative outside i?isco if the s P m of black hole q is sufficiently large, firstly pointed out 
by Tagoshi et al. ll38tl . 

T[0l 

To see how serious the problem is, we revisit the Taylor flux in the test particle limit. The spin independent terms in £ nP N have 
been analytically calculated up to 22PN order and the spin dependent terms up to 8PN order by Fujita lT42l l43ll . The expression 
of the normalized 8PN energy flux including complete spin dependent terms is schematically expanded as 5 

£$Sk*,9) = EE i( "' p) (9)^aogwr , (id 

n=0 p=0 

where p max is the maximum integer that does not exceed n/3. In FIG.[T] we depicted L^ p |,(x, q), the normalized nPN Taylor 
flux in the test particle limit ( [Til truncated at various PN orders for q = 0.9. The curves terminate at the ISCO frequency. We 
also plot the exact numerical energy flux in the test particle limit as a reference. The exact numerical energy flux is manifestly 
positive definite for x < £isco(0.9), while all the nPN Taylor fluxes in FIG. [T]cross zero at some x — Xo(q) < &isco(0-9). 
In addition to the already known 2.5PN and 4PN cases 6 [38], we also find that 5PN, 5.5PN, 6.5PN and 8PN fluxes become 
negative before x reaches a;isco('?)- Even for a moderate value such as q = 0.7, we still observe the flux to cross zero, say at 
8PN order. Once this happens, the integrand of © diverges, and then the phase evaluated by using the truncated flux £^p°]j does 
not make sense. For a reliable extrapolation, it is necessary to use a resummed expression for the energy flux that gives at least 
a finite estimate of the phase for the truncation at any PN order. 



5 The explicit expression for £ will be made available elsewhere in an appropriate form l43ll . 

6 While Tagoshi et al. l38ll reported that 3PN Taylor flux with q = 0.9 became negative for small radius, we find that the 3PN Taylor flux is always positive 
definite for all radius irrespective of the spin parameter of the Kerr black hole. 
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B. The exponential resummation and an improved hybrid energy flux 

To overcome the drawback of the nPN Taylor flux, a naive requirement will be to guarantee the energy flux to be always 
positive by resummation. Among various known resummation techniques JUrH^], to the best of our knowledge, only the 
factorized resummation ensures the positivity of the energy flux. This resummation was proposed by Damour et al. IdTII for the 
equal mass non-spinning binaries in a circular orbit. Fujita and Iyer iBijill and Pan et al. |[59ll applied it to the test particle in a 
circular orbit around a Schwarzschild and a Kerr black hole, respectively. Here, we propose another even simpler resummation, 
which is defined by 

32 

£°p P N (x,q) := yj/ 2 a; 5 exp[i°p P N (a;,(7)] , (12) 



with 



(13) 

truncated at nPN order 

where n' > n is understood. We refer to this resummation and the energy flux ( fT2l as the exponential resummation and "the 
?iPN exponential resummed flux", respectively. This flux is manifestly positive and can incorporate any higher order corrections 
if the counterparts in the nPN Taylor flux is given. Furthermore, it significantly improves the accuracy of approximation l43ll . 
For definiteness, we explicitly write the definition of 'C^n"' an( l ^kPN*' as 

Cpf ] (*,9) f ^Vexp^M)] , Cpf W) := (f ^ "P^Sf W)]) I%g*{x,q) , (14) 

with L° p P N = i^p P N°' + ^£«pn + ■ • • ■ We show the explicit form only for L° n x p p ^ as 



rexp[1] 35 17 3/2 30523 2 / 136229 101 \ 5/2 

j„ P v (x.q) = x H qx ' H x + q ir ]x ' 

nPN V . H) 12 T- 4 H 4032 I 4032 y 



43670915 41 2 tt \ , . 70075 



)x 3 + 7 -^™ 7/2 + 0(x\ q x^,qi). (15) 



12192768 48 ^ 2 9 

Combining the exact numerical flux obtained in the test particle limit, £|"' n , and the PN flux, £pn, one can obtain a better 

[01 T[ll 

estimate for the energy flux. In the case of the Taylor flux, these two are combined simply by summing £f u J n and '. In the 
case of the exponential resummed energy flux (1121) . we need to do this summation at the level of the exponent, L cxp . Namely, 
expanding L^p p N in powers of v, we replace L^}^ 1 with the one corresponding to the numerical flux. Then, we obtain 



£ nPN ( x > l) = 4uii ( x > l) exp vL^ 1] (x, q) 



(16) 



which we call "the nPN hybrid flux". 

In the analogous way as before, we introduce "the known nPN exponential resummed flux" and "the known nPN hybrid flux", 
and distinguish them from their respective counterparts by associating " ~ ". Here the truncation for the known part is made at 
the level of L oxp W . Namely, £^pjl is truncated at the same PN order for each order of spin that is included in the known terms 
^nPN- F° r me exponential resummed flux, we can also introduce £^pn°'' * :rie counterpart of 'C^n 1 ' in the test particle limit, 
as in the case of <£^ P °]!j. Namely, i°p P N ' is truncated at the same PN order for each order of spin that is included in the known 

terms ^pn- On the other hand, i^p 1 ^ ' does not make sense. In the following we identify the hybrid flux with the exponential 
resummed flux for the leading order in the mass ratio. We summarize our notation for various fluxes used in the rest of our paper 
in Table H] for readability. 



IV. THE SCALING LAW OF THE COEFFICIENTS IN THE ENERGY FLUX 

Now we tackle our main issue: how to evaluate the magnitude of the yet unknown part of the energy flux, 5PN - (When 
we discuss general issues independent of the form of the energy flux, we simply suppress the labels "T", "exp" and "hyb".) Our 
strategy here is to establish the scaling property in the coefficients of the PN expansion of the energy flux in the test particle limit 
for fixed q. (See also Ref. l44ll .) The argument of the nPN exponential resummed flux, L° n X pj^\ is expanded as 



Xh(q) 



(17) 



7 



The flux truncated at rtPN order 
The residual part of £|jp N 
The known part of flux truncated at nPN order 



The residual of £ 



The normalized nPN Taylor flux 
The exponent of nPN exponential resummed flux 



Symbol 

fli] 
^nPN 

rH 

fW 
^nPN 

fli] 

rT[i] 

^nPN 
rexp[i] 
^nPN 



TABLE I: Our notation for various fluxes in this paper. Upper index [i] refers to the order of truncation with respect to the mass ratio v. 
L = £ [0] + uL [1] + 0(u 2 ). We also put label "T", ' 'exp", '"hyb" and "full" which respectively correspond to Taylor, exponential, hybrid and 
numerical complete flux, to distinguish the flux type. 
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FIG. 2: The absolute value of the contribution to the energy flux from each PN order |CVi,o| at x = x\x(q), with the horizontal axis being the 
PN order n. Note that the values of |C n ,o| are accidentaly small with some fixed q. We safely exluded these values from the plot. 



where 



2/3 



x h (q) :-- 



p3/ 2 , „ 



(18) 



with 



1 + cos [ — arccos q 



(19) 



is the value of x corresponding to the radius of the circular orbit on the light ring 14811 . 

It is known that for circular orbits the source term of the Teukolsky equation has a simple pole when the orbit is on the light 
ring radius, i.e., at x — x\ r (q) ll60ll6lll . The light ring radius is the innermost radius for the presence of an (unstable) circular 
orbit, and there the particle energy per unit rest mass diverges. As a result, the energy flux also diverges there. Besides the light 
ring singularity, another origin of singularity in solving the Teukolsky equation is the presence of quasi-normal mode poles, 
which are zeros of the Wronskian. However, the absolute value of the first quasi-normal mode frequency, which is the smallest, 
is almost identical to the frequency of GWs from the source at the light ring ll62tl . Hence, we expect that x\ T determines, at least 
approximately, the convergence radius when the energy flux is seen as an analytic function of x. If the PN expansion of the 
energy flux converges as long as x < xi r (q), x = x\ T (q) is the first place where the series ceases to converge. 

With this expectation in mind, we plot the contribution from each PN order with x = x\ T (q), |C n ,o|> in FIG. G] for various 
values of q. Here q = 0.998, is a bound for possible maximum spin of an astrophysical black hole called Thorne limit l63ll . As 



is expected, |C n o| is roughly independent of the PN order n, as typical at the boundary of convergence for \q\ < 0.9. Contrary 
to our expectation, for \q\ > 0.9, C nj o| continues to increase up to 8PN, and the series seems to diverge. However, this does 
not mean that PN series does not converge for x ps x\ T (q) with \q\ ps 1. When we make a similar plot for each partial wave 
contribution for a higher multipole component, the contribution to the energy flux peaks around a higher PN order, but the 
series is convergent, although we cannot check extremely high multipole components. On the other hand, the summation over 
multipole components is known to converge from the numerical computation of the energy flux by Fujita and Tagoshi l5ll l52tl . 
These facts may indicate that the convergence depends on the order to take summation. The point of FIG.[2]is that |C ni o| shows 
a nice scaling property with respect to the PN order for any value of q. Even for q close to the extremal limit, the plot becomes 
flat if we substitute a little smaller value of x instead of x\ T . 

Now we are in the position to discuss how to estimate the order of magnitude of the yet unknown part of the energy flux, 
i">3 5PN - The second order perturbed Einstein equation would schematically take the form Oh^ = (T[z + Sz] — T[z\) + 
h^T[z] + (VM 1 )) 2 , where and are, respectively, the first and the second order perturbations induced by a body. Here 
z is the background Kerr geodesic and Sz is the 0{v) correction to it l65ll . We know that all the first order perturbations have 
singularity only at x « X\ T (q). The Green's function is basically the same as in the linear case. Hence, even for the second 
order perturbation, the convergence radius in x will be the same as in the linear case, and thus a similar scaling for the PN 
coefficients will be expected for the next leading order in the mass ratio, too. Then, we can guess the amplitude of the higher PN 
coefficients from the first few terms in the PN expansion that are known from the standard PN calculation. 

One may worry that has singularities at the quasi-normal mode frequencies. When we consider the second order per- 
turbation, there may arise high frequency components which are absent at the level of linear perturbation. However, this might 
not happen. The point is that the metric perturbation caused by a circular geodesic has a helical Killing vector 1 . At the level 
of second order source term, this helical symmetry is broken due to the presence of the deviation from the geodesic, Sz. Apart 
from this contribution, however, the source term keeps the helical symmetry and hence the frequency that appears in a partial 
wave labelled by (£, m) is mil only. Namely, higher frequency modes do not arise. Even if we take into account Sz, the time 
scale associated with Sz is as slow as il^ 1 x 0(i , ~ 1 ). Therefore, the presence of Sz alters the frequency mil only by a small 
amount of 0(vil). This essentially does not change the convergence radius in x. 

Even if we get the order of magnitude of the PN coefficients, the scaling property tells us nothing about the actual sign of each 
term. Therefore, when x is close to x\ T {q), it becomes difficult to guess the order of the magnitude of the infinite summation. 
However, since the higher order PN terms are suppressed by the power of x/x\ r (q) , the summation is dominated by a few leading 
terms when xJx\ T {q) is reasonably small. In that case, we do not have to worry about the infinite summation. Here, we bravely 

step a little forward by making the following proposal that the ratio between £^p N and £>'„ PN will be the same order as the 
ratio between £^ PN — ^ora anc ' ^>Lpn- B ase d on this assumption, we extrapolate the known results to estimate the unknown 

^>nPN' 

V. THE DEPHASING DUE TO THE AVERAGED DISSIPATIVE PORTION OF THE SECOND ORDER SELF-FORCES 

The main focus of this section is to evaluate the dephasing from the averaged dissipative part of the second order self-forces 
in E(I)MRJs, based on the idea of the extrapolation proposed in the preceding section. In this paper, we measure the magnitude 
of the dephasing due to higher order corrections to the energy flux SL by 

A^[L,SL}:=-2 dx - * W 

Here the absolute value of the difference of fluxes is taken in the integrand. The factor (£ + (5£ ) ~ 1 — £ ~ 1 can change its signature 
in the domain of the integral for some parameter region of E(I)MRIs. If we do not take the absolute value, therefore, there might 
be an accidental cancellation between positive and negative contributions. Then, the order of magnitude of the evaluated A$ can 
largely deviate from what we really want to measure. Notice that, even if A$ defined without taking the absolute value in the 
integrand strictly vanishes, the deviation in the GW waveform is still detectable. Therefore, to avoid this possible underestimate 
of the dephasing, we take the absolute value of the difference. 

As a general remark, we would like to mention the mass dependences of < 3>[<cj°f], the phase for a certain reference flux 
which includes the leading PN terms in the test particle limit, and A$[,c[.^ f , SL^] and A$[£^ f , <5£M], the dephasings due to 



L{x)+SL{x) L{x) 



(20) 



7 The helical Killing vector t a + f2</> Q , that remains to generate a symmetry even for the perturbed spacetime: where t a and (f> a are the asymptotically time 
translation and rotational Killing vectors of the background Kerr black hole, respectively. See e.g.Ref. ^6$\ for more detailes. 
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FIG. 3: The dephasing due to the higher PN corrections in the test particle limit, i/A$[£™jfj!j , ^^pj,] as a function of v/M for q = ±0.9. 

5L I°l and 5L M , which are some portions of the energy flux at the leading order and the next leading order in the mass ratio, 
respectively. From the integral in Eq. (|4j, which defines xq, one can factor out M-qyi/v since 

oc AW, ^oci/ 8 , (21) 

for a given x. Hence, we find that xq is a function of Msr/v. Substituting 5L<- ' oc v 2 and SL^ oc v 3 together with the 

relations (ED into the definitions of v X v x A$[£[.° ] f , SL^] and A$[£™ , SL^}, we find that their mass dependences 

remain only through xq. Therefore, one can conclude that v x *[«C™f]' v x A$[£j° ] f , SL^} and A3>[£[°J f , £W] depend on Af BH 
and only through the combination Msa/v for a large mass ratio i/ < 1 for a fixed observation period before the plunge. 

A. The error caused by the post-Newtonian truncation in the test particle limit 

Before we evaluate the dephasing due to the yet unknown part of the energy flux at the next leading order in the mass ratio, 
we will assess the magnitude of the dephasing due to the PN truncation in the test particle limit A$[£^*p N , £>J lPN ] in this 
subsection and that coming from the known PN terms, A$[£|"jj, £^p N ], in the succeeding subsection. 

Here, we study the quantity A$[£j i ° pN , £^| iPN ] to examine the PN convergence in terms of the phase error for various 
E(I)MRIs. We show the results only for the exponential resummed flux because the phase for the PN Taylor flux becomes 
ill-defined for several PN orders. (The hybrid flux is identical to the exponential resummed one by definition in the test particle 
limit.) 

In FIG. |3] we plot vA^[L^^\ ^™^pn] as a function of v/M for q = ±0.9. The main message of these plots is that the 
PN convergence is almost uniform except for the case of large black hole spin for large v/M. When q is large, ISCO becomes 
close to the light ring radius where PN convergence becomes very slow. Even in that case, for lower PN orders the convergence 
is rather smooth for large v/M. This is because up to 2PN order the phase correction is dominated by the lower bound of the 
x-integral in d20l i. if the initial separation of the binary at the time 1 year before the plunge is sufficiently large. The reason why 
we focus on large v/M is, as will be discussed in Append. M the dephasing from the known PN terms is already well suppressed 
below 1 rad if v/M is small enough. Thus, the PN convergence for small v/M does not affect the following discussion. 

Notice also that xq becomes smaller and smaller for larger v/M. The 2.5PN case is marginal in the sense that all range of x 
almost equally contributes. By contrast, the contribution near the ISCO dominates for the corrections at the 3PN order or higher. 
Although the PN convergence is not clearly seen for q > 0.9, we should note that the phase error is never extremely enhanced 
for some particular post-Newtonian order, largely exceeding the values for 3PN or 3.5PN order, which we mainly focus on in 
the following discussion. 

B. The expected dephasing from the unkown higher order post-Newtonian terms 

Now we move on to our main issue: evaluating the dephasing due to the yet unknown higher order PN terms, A$[£|^ n ± 
^3 5PN' ^>3 spn]> based on the idea for the extrapolation proposed in Sec. [TV] The scaling argument will tell us that the dephas- 
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FIG. 4: The expected residual dephasing A$ gue sa [-Cfun + -^[Jpni ^>Lpn] caused by the unknown part of the averaged dissipative second 
order self-forces. For n — 3, the 3.5PN known spin independent part of the flux at the next leading order in the mass ratio is treated as unkown. 



ings coming from the respective PN terms at the leading order in the mass ratio will be roughly proportional to the corresponding 
next leading order dephasings. Namely, 

Aar/M fM rM l • A*[f [0] r [0] i~A*r? [o1 r [1] i-A<j>rr [o1 +r [1] r [1] l c?7\ 

^^^OPAf'^nPN ~~ ^OPNl ■ ^^L^nPN' ^»iPnJ ~ ^^L^fuUi ^nPNJ ■ ^^l^full + ^nPN' ^>nPNJ i \ LL > 

will hold. Recall that £^ PN consists of the terms at the leading order in mass ratio up to the same nPN order that is available in 

'WiPN- 

Here, we have basically three options in the choice of the flux among the Taylor flux, the exponential resummed flux and 
the hybrid flux. (Recall that the hybrid flux is defined as such that is identical to the exponential resummed flux at the leading 
order in the mass ratio.) On one hand, the Taylor flux is problematic since it can be negative before reaching ISCO for some PN 
orders, say 2.5PN, with modelate value of the spin parameter as we mentioned earlier. In such cases the phase before the plunge 
is not well defined. On the other hand, the difference between the exponential resummed flux and the hybrid flux is negligible. 
Therefore, we here consider the hybrid flux only. 

The residual dephasing estimated by using d22b 

A* \r [0] -4-/-W r [1] l— ^^»PN'^>nPNJ AfrifM r [1] 1 nv> 

with n = 3 and n = 3.5 is depicted in FIG.IU The trend of the estimated dephasing is independent of the truncated PN order. 
The point of FIG.|4]is that the residual dephasing is rather supressed over the whole range of the bainary parameters. Indeed, the 
residual dephasing is at most about 10 rad in the case of 3PN truncation. In the case of 3.5PN truncation, the maximum value of 
the residual dephasing is bigger by factor 2 or so. Since the plot is given in the unit of radian, the value must be divided by 27T to 
translate it into the number of cycles. From the above results, one may say that the residual dephasing due to the yet unknown 

PN corrections at the next leading order in the mass ratio, A$ gucss [,c[°' n + £^ PN , £>Lpn] ls not negligible, adopting one cycle 
as the criteria for the significant dephasing. 

In FIG.|U we see a tiny bump around vMq/M m 10~ 10 — 10~ n for some values of q. The reason why this bump appears 

at this position can be understood from the observation that £|^ PN — ^opn crosses zero within the domain of the integral, 

xq < x < xisco when the bump appearers. When £|°jl N — ^opn crosses zero, the factor in the denominator of Eq 

nPN ~ ^opjvI suffers from more or less accidental suppression even though we take the absolute value i 

integrand of A$ in (l20b . We denote the value of x at which the flux £^ PN — £[ip N crosses zero by a; cross . Then, if xq <C x 



cross 



or xq 3> a; cross , this suppression does not produce much effect on the estimate of A$[Lqp N , £„ P n — ^opjv]- Therefore, the 
suppression becomes significant only for xq rj across- The value of x cross is rather close to the value at ISCO, £isco> but 
typically not extremely close to it. Hence, xo « x cross occurs when xq is neither extremely close to ISCO nor very small like 
xq <C 1, which corresponds to the parameter region vMq /M w 10~ 10 — 10~ u . Since now we find that this tiny bump is to be 

attributed to an accidental zero in £j°p N — £npjv> tn is bump would be regarded as an artificial feature. If it is fare to remove the 
bumps from FIG.|4j we will find that the significantly large dephasing (in the sense of exceeding 1 cycle) will be expected only 
for vMq/M > 10" 9 . 
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FIG. 5: The expected residual dephasing A<&g ue ss[<Cf u n + ^Lpni ^>1iPn] wnen we assume that the spin dependent terms at the lower PN 
orders are all known. 



C. Unknown dephasing expected when we know the lower PN order non-linear spin dependent terms in the energy flux 

In the preceding subsection, we find that the dephasing caused by the unknown PN higher order terms at the next leading 
order in mass ratio can be 0(10) rad or more. However, these unkown terms in the flux include non-linear spin terms at 3PN 

or lower order and all spin dependent terms at 3.5PN order, and hence we excluded the corresponding terms from £[[p N , which 
appears in the expression d23b . In FIG.2J relatively large dephasing is observed for large for which the non-linear spin terms 
will be important. 

Here we consider how the estimate presented in the preceding subsection is modified once we obtain all the spin dependent 
terms up to 3PN or 3.5PN order. What we need to evaluate is the expression obtained by removing " ~ " from d23l , i.e., 

\ 3,\ AO] A°] 1 

a* \r l0] +r [1] £ [1] i— Aq ^»PN>^>nPNl A^rr [0] r [1] i r?/n 

Here, one may think that £^p N in the factor A$[£^p N , £jj PN ] should have been also replaced with £„pN' Since we do not 

have the expression for £^p N at hand, we cannot perform this replacement. However, the difference between £^p N and £Pp N 
will not be significant because their dominant PN terms are common. 

The resultant dephasing A$ gucss [£|^ n + £|Jp N , £>Lpn1 i s depicted in FIG. [5] In contrast to the FIG. 2] except for the near 
extremal spin case (q = 0.998) with the 3.5PN truncation (n = 3.5), the residual dephasing for large v/M evaluated by using 
the hybrid flux does not exceed 1 rad. The residual dephasing exceeds 1 rad only at around the peak, which is likely to be an 

artifact due to the suppression in A$>[L^p N , <&|° PN — ^opjvl similarly to the previous case. The above result indicates that the 
dominant source of the error caused by the unknown part of the flux is in the spin dependent terms at the lower PN orders up 
to 3PN or 3.5PN order. Namely, the knowledge of non-linear spin dependent terms in the energy flux is crucial for reducing 
the uncontrolled dephasing caused by the averaged dissipative second order self-forces. In addition, once we obtain these spin 
dependent terms of the energy flux, the exponential resummation will improve the accuracy of it to the level almost sufficient 
for detection of E(I)MRJs in the whole interesting parameter region. 



VI. SUMMARY AND CONCLUSION 



To evaluate the second order self-forces in the context of the black hole perturbation is a challenging issue motivated by the 
goal to extract information from extreme (intermediate) mass ratio inspirals (E(I)MRIs). In this paper we have assessed the 
dephasing of the GW waveform from quasi-circular E(I)MRIs, caused by the averaged dissipative part of the second order self- 
forces, focusing on giving an order of magnitude estimate of the influence of the yet unknown higher order post Newtonian (PN) 
corrections. Using the balance argument, the dephasing is related to the correction to the the emitted energy flux. Although it will 
still require further efforts to calculate the second order self-forces using the black hole perturbation, they are partially already 
known from the standard PN calculation. We gave an estimate on how large dephasing will be caused by the yet unknown higher 
order PN terms for the last 1 year of inspiral before the plunge, exhaustively exploring the whole possible parameter region of 
E(I)MRIs. 
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To give a guess for the yet unknown higher order PN terms, we first introduced a simple resummation method for the energy 
flux, which we call the exponential resummation. This is simply obtained by exponentiating the energy flux and truncating it at 
the known PN order in the exponent. This resummation has three merits: it ensures the positivity of the energy flux, it accelerates 
the PN convergence, and it can be applicable if we just know the PN Taylor flux. Since the PN Taylor flux in the test particle 
limit with large q is known to be negative outside the innermost stable circular orbit (ISCO) at several PN orders, the total phase 
before the plunge becomes ill-defined. Hence, we cannot discuss the amplitude of the higher PN corrections by using the PN 
Taylor flux. The exponential resummation solves this problem and improves the PN convergence in terms of the total phase 
before the plunge. When we discuss the finite mass corrections to the energy flux, we can combine the idea of the exponential 
resummed flux with the known exact energy flux in the test particle limit. We also proposed such a phenomenological energy 
flux, which we call the hybrid flux. 

To examine the dephasing from the unknown part of the averaged dissipative second order self-forces, we need some extrapo- 
lation. Using the brand new 8PN energy flux in the test particle limit [43] that includes all spin dependent terms, we discovered 
that the order of the magnitude of the absolute values of the coefficients up to 8PN order approximately follows a simple scaling 
law, which is what we expect when the actual convergence radius of the PN expansion is at the light ring radius, irrespective of 
the value of dimensionless spin parameter q. Since there is no reason to expect that the energy flux diverges outside the light 
ring radius even if we take into account the finite mass corrections for quasi-circular E(I)MRIs, we are motivated to assume that 
the ratio between the magnitude of the terms at 0(i>°) and O(v) at the same PN order will roughly stay independent of the PN 
order. 

Based on this assumption, we estimated the unknown portion of the energy flux that comes from the higher PN terms at the 
next leading order in the mass ratio and evaluated the residual dephasing due to them. We find that the residual dephasing may 
exceed one cycle for vMq/M > 10~ 10 and the spin of the Kerr blackhole \q\ > 0.5, assuming 1 year observation period. 
For some parameters, we found a little enhancement of the estimated dephasing at vMq/M » 10 -10 — 10 -11 . Since this 
enhancement is likely to be attributed to a mild accidental cancellation in the factor in the denominator of the estimator of the 
residual dephasing that we adopted, it might be fair to remove the bump that arises for this reason. Then, the residual dephasing 
exceeds one cycle only for relatively large 0(v/M) with \q\ > 0.5. Even for rather extreme case like q = 0.998 (Thorne limit), 
the expected residual dephasing is at most a few cycles or so. 

In the unknown flux, non-linear spin terms at lower PN orders are also included since they are not yet calculated in the context 
of PN approximation even at the leading order in the PN expansion. As large residual dephasing is expected only for large \q\, 
one may suspect that the residual dephasing might be dramatically reduced once we know all the non-linear spin dependent terms 
up to 3PN or 3.5PN order. Therefore we made analogous plots for the expected residual dephasing assuming that the non-linear 
spin dependent terms are also available. We found that the expected residual dephasing is further suppressed to be less than 1 
rad for \q\ < 0.9, as we expected. Even in the extreme case with q = 0.998, the expected dephasing becomes less than 1 rad if 
we use 3PN truncation to evaluate the residual. We think the reason why 3.5PN truncation is worse is simply because, for large q 
close to unity, ISCO is very close to the light ring radius and hence the PN convergence is very poor. Thus, the resulting estimate 
of the residual also fluctuates by a large order of magnitude in this limit. However, even if we rely on the 3.5PN truncation, the 
expected dephasing is at most about 1 cycle for q = 0.998. 

Although one year observation period is assumed in all the above estimates, the results do not depend much on it. The 
change in the observation time varies the relation between xq and v/M, and hence the curves in FIGs.[4]and[5]are shifted just 
horizontally leftward. The expected dephasing becomes large only for large v/M, but xq is already sufficiently small in this 
case. Therefore the effect of changing the observation period for large v/M appears only in FIG. 2J in which lower order PN 
terms are contributing. In fact, one can see that the curves in FIGs.[5]are already saturated for large v/M . 

To conclude, the residual dephasing caused by the unknown averaged dissipative part of the second order self-forces is 
estimated to be at most a few cycles. This dephasing is mostly to be attributed to the lower PN order non-linear spin dependent 
terms. If the PN expansion to the 3PN or 3.5PN order is completed including the dependence on the black hole spin, the residual 
is expected to be further reduced. 

As a final remark, we should note that there is a possibility that we might be underestimating the residual dephasing by an 
order of magnitude, since the PN convergence is not very smooth, especially for positive large q. Even though our estimate turns 
out to be a good approximation to the real magnitude, we still have a chance to have a golden event with a large signal to noise 
ratio. In that case, one cycle or even one radian may not be a sufficient accuracy that is required for the template waveforms. 
Then, the dissipative part of the second order self-forces that cannot be captured by the PN expansion becomes necessary to 
extract the best physics from the observational data. Also we should mention that we focused on quasi-circular E(I)MRIs in this 
paper, but our scaling argument will not apply any more for significantly eccentric orbits, since there are many modes whose 
frequency exceeds the naive PN convergence radius that will be given by the absolute value of the first complex quasi-normal 
mode frequency. Therefore, the hybrid energy flux cannot be expected to remain a good approximation near the plunge for such 
systems. Hence, our analysis does not at all discourage the study on the averaged dissipative part of the second order self-forces. 
Our claim is that it will be possible to perform nearly the best analyses of most of (nearly) quasi-circular E(I)MRIs, without 
waiting for the full development of our knowledge about the averaged dissipative part of the second order self-forces. 
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FIG. 6: The dephasing caused by the known terms up to the 3PN order at the next leading order in the mass ratio, A$pn[<c[°ii, "C^pn ]> f° r 
various dimensionless spin parameters q of a Kerr black hole, based on the 3PN hybrid flux U6I >, 
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Appendix A: Dephasing from the known post-Newtonian terms 

This Appendix is dedicated to see the dephasing from the available PN terms in the next leading order in the mass ratio. In the 
bulk of our paper, we introduced three different types of PN flux formulas, Taylor, exponential resummed and hybrid, but there is 

no qualitative difference among A$ PN [-C|°n, £^pnL A<f> PN [,£|°' n , ^ P n'] an ^ A^pnI^Iuh, ^nPN^l' as l° n § as we c h° ose me 
exact numerical flux <cM, as the reference flux. This is because the significant difference among various fluxes appears only in 

the higher PN residual parts such that we discussed in the Sec. [V] Therefore, we only show A$pn[£[°j[, ^pn^I m FIG- El We 
stress that what we show is the dephasing caused by the terms not at the 3PN order but up to the 3PN order. The main message 
of FIG. |6] is that the corrections to the phase of GWs due to the terms higher order in the mass ratio are suppressed below 1 
rad, even if we include the leading order of the PN expansion, for a binary with vM & /M < 10~ 12 , say, when the mass of the 
central Kerr black hole is greater than 3 x 10 6 M Q with the satellite mass fixed to \QM®- This holds irrespective of the value of 
the dimensionless spin parameter q. This plot clearly denies the naive statement: the effects of higher PN corrections should be 
important for a more massive central black hole because the satellite stays near the ISCO for a longer period. 

As for the results in literature to compare, for the binary with A/bh := lO 6 A/0, q = and v = 10~ 5 , Heurte and Gair ll35ll 
reported 1.5 rad dephasing after the last one year of inspiral based on the 2PN Taylor flux, while Yunes et al. Ik36tl reported 0.71 
rad dephasing after the last 1 1.5 months of inspiral for the same binary parameters. Following our definition of the dephasing, 

we calculate A<&pn[<c[uh, ^pn 1 '] t0 ^ n< ^ ^ 1-46 rad, which we think is in a good agreement with the previous results within the 
variance due to the different definition of dephasing. 
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FIG. 7: The dephasing ^A<E>[£f ull + L 2 5PN , spn! due t0 tne horizon absorption flux and the dephasing i/A$[£°^ , .-C^pn] due to 
the higher PN corrections, for n = 3 and n = 3.5, and q = ±0.9 



Appendix B: Dephasing due to the horizon absorption flux 

In the bulk of our manuscript, we neglected the flux absorbed by the horizon. In the test particle limit, Tagoshi et al. |6f 
and Yunes et al. l36ll have already shown that the energy flux absorbed through the horizon cannot be neglected, leading to a 
large dephasing, especially when the Kerr black hole has large spin. However, this does not mean that the corrections due to 
the horizon absorption flux at the next leading order in the mass ratio are also important. Naively the standard PN formalism 
is not suitable for calculating the horizon absorption flux since the black hole horizon is beyond the reach of the standard PN 
expansion, alt houg h there is a direction to evaluate the horizon absorption flux by relating it with the tidal field around each 
black hole 167116*1 , 

In this appendix, we shortly address the dephasing caused by the flux absorbed through the horizon in the test particle limit 
to see what can be said on the dephasing due to the horizon absorption flux at the next leading order in the mass ratio. The flux 
absorbed through the horizon in the test particle limit is expanded as l66ll 

4°p H n ] W - f^ 5 {(-^-^ 3 )^ 5/2 +(-g-|9 3 )- 7/2 + o(x 4 )} , (Bi) 

and the exact numerical flux is also calculable OH]. We denote the contribution from the first term in the curly brackets in 

the above expression by ^^spn an d me remaining part of the horizon flux by -C^^spN' m a simile manner as before. As we 
mentioned above, the horizon absorption flux at the next leading order in the mass ratio is yet under development, but the part that 
requires the second order metric perturbation should start with 1PN order higher than the leading terms in the test particle limit, 
starting at 3.5PN order. Hence, we adopt L > 2 H 5PN' excluding £ 2 ° 5 p' N , as tne counterpart of the unknown horizon absorption flux 
at the next leading order in the mass ratio. Under this consideration, we compare the dephasing due to the horizon abosorption 
flux 

^^I'Wull + ^2.5PN' ^>2.5PnJ y aA > 

to that coming from the PN truncation A$[£°p P f[f ' , -C^p ^] for n — 3 or n = 3.5 in the test partilcle limit. The results are 
plotted in FIG.|7]for q = ±0.9. 

FIG.|7]indicates that j/A$[£[°jj ± ^^sra, ^>2 spn! stavs at most O(10 _1 ) rad for all E(I)MRJ parameter region. Indeed, it 
is, at most, about 1/10 of i^Af&^^p 1 ^ ' , -C^p ^]- F° r large q, the smallness of the effect of the horizon flux relative to the higher 
order PN corrections can be understood from the work by Hughes 16^1 . who numerically found that the energy flux absorbed by 
a Kerr black hole is at best 10 times smaller than the one emitted to the infinity, with the aid of the fact that the PN expansion 
shows a poor convergence near the ISCO for large q. Thus, we expect that the residual dephasing due to the horizon absorption 
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flux will be minor compared with the higher PN flux to the infinity, also in the next leading order in the mass ratio. 
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